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Abstract 

We extend the exactly solvable Hamiltonian describing / quantum oscilla- 
tors considered recently by J. Dorignac et al. by means of a new interaction 
which we choose as quasi exactly solvable. The properties of the spectrum of 
this new Hamiltonian are studied as function of the new coupling constant. 
This Hamiltonian as well as the original one are also related to adequate Lie 
structures. 
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1 Introduction 



A few years ago the Bose-Hubbard (BH) model describing a non linear optic model 
with / bosonic oscillators in interaction was studied in some details [1],[2],[3]. In 
particular it was revealed that the spectrum of the restriction of the Hamiltonian to 
the subspace of vectors involving two quanta has a remarkable property when the limit 
of large / is considered : the spectrum separates into two pieces, one line of discrete 
states forming the so called "sohton band" and another region where the eigenvalues 
form a continuum. The total Hamiltonian describing the BH model commutes with 
the operator counting the number of quanta. As a consequence, the Fock space 
V of quantum states can be decomposed into an infinite flag of flnite-dimensional 
subspaces V^, n = 0, 1, 2, . . . that are left invariant by the Hamiltonian and the whole 
spectrum can be constructed algebraically. Recently, these ideas were generalized to 
exactly solvable Hamiltonians with higher orders in the particle operators and it was 
shown that similar results hold. 

Such properties are intimately connected to specific features of Lie structures. 
More precisely, we prove in Section 2 that the BH model is the sum of scaling opera- 
tors of the Lie algebra sl{f; R) supplemented by bihnear fonctions of the (diagonal) 
operators generating the Cartan subalgebra of this structure. Thus, the number / of 
sites fixes the Lie structure subtended by the model. We also prove the fundamental 
role played by the number n of quanta picking up the irrreducible representations 
(irreps) that are concerned with. 

Operators enjoying the algebraic property above are called exactly solvable [4]. 
The family of exactly solvable Hamiltonians is rather small; however if we replace 
the requirement that the Hamiltonian preserves an infinite fiag of finite-dimensional 
subspaces by the weaker requirement that one finite dimensional subspace is preserved 
by the Hamiltonian, we are left with the notion of quasi-exactly-solvable (QES) [5] 
operators and/or equations. In this case not all but a finite part of the spectrum can 
be computed algebraically. 

In Section 3, we apply the ideas of QES equations and we extend the exactly 
solvable BH Hamiltonian by a new term which preserves only a finite-dimensional 
subspace of the Fock space V containing the subspace V2 where the splitting of the 
eigenvalues occurs in the normal BH model. The new term is characterized by a new 
coupling constant, say A. 

Then, in Section 4, we relate the operators involved in this new model to generators 
of an ad-hoc Lie structure namely the Lie orthosymplectic superalgebra osp{l/2f; R). 

Section 5 is devoted to the spectrum of this QES operator in function of A. The 
analysis of the eigenvalues can be achieved along the same lines as in [1] and our 
results suggest that the splitting between the soliton band and the continuum still 
occurs for the QES Hamitonian. We illustrate these results by some examples in 
Section 6. 



2 



2 Group theoretical approach of the BH model 



The BH Hamiltonian is [1, 2] 

/ 7 
Hbh — — '^[(^jO.j+i + OjOj-i + T^aJaJajOj] (1) 

where the bosonic lowering and raising operators Uj, Oj obey the usual commutation 
rules [ ]] = 5ij and the following periodic conditions : 

a/+i = oi , a\^i = al- (2) 

The main property of the operator Hbh is that it preserves separately any subspace 
Vn of the Hilbert space with n quanta, i.e. the space generated by the vectors of the 
form \ni, 712, ... , n/) for fixed n = rii + n2 + ■ ■ ■ + rif. 

Let us now consider the group theoretical approach of such a model. To do so, 
we fix the grading of the operator Ojak as (j — k). The interacting part of the BH 
Hamiltonian is thus a linear combination of (/ — 1) operators of grading +1 and 
(/ — 1) other operators of grading —1. Following the usual Lie bracket as well as the 
grading rule 

[grading{j), grading{k)] — grading{j + k), (3) 

the operators of grading +1 i.e. a^aj_i give rise to (/ — 2) operators of grading +2 
namely ajaj_2 while the operators of grading —1 i.e. 0]%+! generate (/ — 2) operators 
of grading —2 given by The process goes on and on until the operators aj-ai 

and a\af of respective gradings (/ — 1) and (1 — /) are reached. To these scaling 
operators, we add the (/ — 1) diagonal (or of grading 0) ones a]^;^aj+i — Ojaj with 

j = 1, ...,/ — 1 and we finally obtain 2{J2jZi i)+/~l = /^~l operators generating 
the Lie algebra sl{f; R). 

Let us now turn to some specific examples. 

The case / = 2 is the first significant one. The Lie algebra sl{2;R) is the one 
subtended by the BH model. It is characterized by the commutation relations 

[Jo, J±] = ±2 J± , [J+, J_] = Jo (4) 

and the three generators are realized through 

Jo = 0202 — a\ai , J+ = a^di , J- = 0102- (5) 

More precisely, a rapid look at the Casimir operator 

C = J+ J_ + i Jo^ - 1 Jo (6) 

can convince us that ^ 

C\ni, n2 >= -n{n + 2) \ni, ^2 > (7) 
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or in other words that the so-called D^^'' irrep of sl{2;R) is under consideration in 
the BH model. 

Such a representation has the feature that it can be realized through the following 
differential operators 

T ^ d T 2 d ^ d 

Jo = 2^- n , J+ — —X - — \-nx , J- — — (8) 

dx dx dx 

while the basis state \n—n2, n2 > is assimilated to the monomial x"'^ forn2 = 0, 1, n. 
The BH Hamiltonian (1) is then 

Hbh = -J+- J- + |(27V-Ar2_ j2) (9) 

= —ix^-f-^ + [{x'^ ~ 1) + ~ ■'■)'^] J nx — ■^n(n — 1). (10) 

where N is the number operator of eigenvalue n. It can be put on a Schrodinger form 
if the change of variables 

x = exp(±^y) (11) 
as well as the " Gauge transformation" 

Tl 1 

X = exp (i-^y) exp (- cosh ^7^)^ (12) 
2 7 

are performed. The resulting potential writes 

1 X T?/ 7?/ 

V{y) = -cosh^^y- (n + l)coshv^y 7-(- - 1). (13) 

7 722 

It coincides with the one studied in [6] (cf. Eq. (170)) and wc refer to this paper for 
the determination of the corresponding eigenvalues and eigenf unctions. 

An operator playing a fundamental role in this knowledge of the eigenvalues and 
eigenfunctions is the translation operator T [1, 2]. It is defined by the property Taj = 
a]_,_iT so that T|ni, 77,2, . . . , n/ >= rii, 712, ... , > . The BH Hamiltonian as 
constructed in Eq. (1) commutes with T . 

Let us consider this crucial operator in the context / = 2. A rapid look at 
Eq. (9) can convince us that in order to commute with the interacting part of the 
BH Hamiltonian, T has to be a function of the sum of the two scahng operators of 
sZ(2; i?) i.e. ( J+ + J_). This function can be specified by asking for the commutation 
of the diagonal part Jq with it on the basis. For instance, if n = 3, the function 
c(J+ + J-)^ — 7c(J+ + .] ) does commute with Hbh when the irrep D'-z) is under 
consideration. The constant c is then fixed according to f.i. the requirement T|3, >= 
|0, 3 > which gives c = |. We thus have 

T = + - 7(J+ + 
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according to Eq. (8). This result is generalized to 
for an arbitrary n. 

Let us now consider the case f — 3. The algebra sl{3; R) is the one subtended 
by the corresponding BH model. It is generated by 8 operators given in terms of 
annihilation and creation operators following the original model or, equivalently, in 
terms of differential operators : 

0201 ~ , 0102 ~ —xldxi—xiX2dx2+nxi , O2O2 — ojoi ~ —2xidxi—X2dx2+n, (16) 

0303 ~ —9x2 ) (^1^2 ~ a;29j:2 +^i^29a;i — nx2 , 0303 — 0302 ~ Xidxj^+2x2dx.2 — n, (17) 

0I5O1 ~ -X2dxj , a\as ~ -xidx^. (18) 

Both forms are two different realizations of the same irrep D{0, n) of s/(3; R) as clear 
from the dimension ^(n + l)(n + 2) of the basis these operators act on as well as the 
eigenvalues of the "hypercharge" and "T3" operators [7] i.e. 

^ = ^(24«3 - - 4«2) ^ X2dx2 - I (19) 

and ^ ^ 

T3 = -(0^02 - ojoi) ~ -xidx^ - -X2dx^ + -, (20) 

respectively. The BH model can thus be written as in Eq. (1) or with the differential 
realization 

Hbh = -j{xldl^ + xldl^ + xiX2dxidx2) + {xl - 1 - X1X2 + j{n - l)xi + X2)dxi 
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jyn — 1 )X2 — x\)Ux2 ~ 't-xi -f 'ILX2 — 

The equivalence of the bases is 



{x\ — 1 — X1X2 — ^{n — l)x2 — Xi)dx2 — nxi + nx2 — -^^{n — 1). (21) 



\n- k - j,k,j >r^ Xl ^ ^xi, j = 0, 1, - A;, A; = 0, 1, (22) 

The embedding of the cases f — 2 and / = 3 is thus clear. It is not possible here 
to convert the BH Hamiltonian into a Schrodinger form but it is rather straightfor- 
ward to determine its eigenvalues and eigenfunctions either on the bosonic or on the 
differential forms. 

Here also, it is possible to express the operator T in terms of differential operators. 
However, it has to be done by hand and no general expression is available. For 
instance, if n = 1 we have 

T=-{xi + x2){dx,+dx2) + l (23) 
or in terms of bosonic operators 

T = o^oi + 0^03 + o|o2 (24) 
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while if n = 2, we obtain 

T = ^(1 + x^X2fdl^ + ^(xi - xlfdl^ + (1 + XiX2){xl - xi)d^,d^, 

-X2{1 + xiX2)dxi + X2{xi - xl)dx2 + xl (25) 

which is not a polynomial of the operator (23). 

The generalization to the cases of higher / does not present any difficulty except 
for its heaviness, technically speaking. 

3 The QES model 

The Hamiltonian which we consider now is given by 

H^Hbh + H^ (26) 
where the second piece is chosen according to 

/ 

Hx = \J2 («1 - 2) + (iV - 2)a,) , (27) 

with the total particle number operator N 

f 

N^Y^A^J- (28) 

If the operator Hbh preserves separately any subspace V„ of the Hilbert space , the 
piece Hx preserves only the subspace Vq ® Vi ® V2 . In this sense, the full Hamiltonian 
H is said quasi exactly solvable [5] since it preserves a finite-dimensional subspace of 
the full Hilbert space. Accordingly, the Hamiltonian can be diagonalized algebraically 
on the subspace. We call this operator the QES BH Hamiltonian. Note that the 
occupation number operator N does not commute with H while T does. 

4 Group theoretical approach of the QES model 

The QES BH model (26) implies that we add now the bosonic annihilation and 
creation operators to their bilinear products previously considered. If we ask for 
commutation relations only, the algebra will not close, that is why we are naturally 
led to associate a Z2-grading to the operators involved in Eq. (26). Hence the 
operators a^- and % have an odd parity (and thus obey anticommutation relations) 
while their bilinear products are even (and have to satisfy commutation relations). 
So, to the previous {p — 1) even operators ajafc, we add now \f{f + 1) even ones 
given by a^aj., their conjugates as well as the number operator which is not a invariant 
anymore. We thus obtain {2p + /) even operators. We complete the structure with 
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the 2/ annihilation and creation odd operators. These operators generate the Lie 
orthosymplectic superalgebra osp{l/2f; R) as is well known [8]. 

The dimension of the osp(l/2/; R) irreps involved by the QES BH model can be 
determined according to the counting of the states : 1 for Vo, / for Vi and |/(/ + 1) 
for V2, giving a total of |(/ + 1)(/ + 2). The situation is thus different from the 
one encountered in the usual BH model. In the BH model, the number / of sites 
determines the Lie algebra while the number n of quanta fixes the involved irrep. In 
the QES version, the number of quanta is fixed at the start and the number of sites 
selects both the Lie superalgebra and its irrep. 

Let us once again put our attention on two significant cases. 

If / = 1, the fundamental irrep, with — as the eigenvalue of the Casimir operator 
[8], of osp{l/2; R) is under consideration. It can be realized through bosonic operators 
as in Eq. (26) or via matricial differential operators. In what concerns the odd 
operators, we have [9] 



a = V^(; = ""V") (29) 



which implies 




[ V2X{2x{x + l)£-x-2 -7(2x^^ + x£) - Ax± 

The equivalence at the level of the states is 

io>^f M,ii>~f M,i2>^( : 1. (31) 




The energies are determined straightforwardly whatever the form of H is. They are 
the solutions of the equation 

+ (7 + 6)^2 + (27 + 8 - Q\^)E - IGA^ - 4A^7 = 0. (32) 

The translation operator T still commutes with H since it commutes with both a and 

. It simply reduces to the identity operator. 
Let us now turn to the / = 2 context. 

The six- dimensional osp(l/4;i?) irrep is under consideration here. It can be real- 
ized with differential operators of two variables [10]. Nevertheless, we choose here a 
differential realization in terms of one variable only since we want to compare the A 
contribution with respect to the original model. In this aim, we consider the Hamil- 
tonian (10) restricted (due to the QES features) to the cases n = 0, 1, 2. It reads 



H 



BH 



/ \ 

2{x'^ -l)^-2x 
Vo ' 2(x2-l)^-4x-7(x^^-x^ + l) y 



(33) 
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At the level of the states we have 



|0,0>^ 







/ 


(")■ 


1,0 >~ 


1 










( " 




,10,2 


>- 


) 









|0,1>^ 




(34) 



Knowing these states as well as the action of the annihilators and the creators on 
them f.i. 



ai|0, >= , ai\l, >= |0, > , ai|0, 1 >= , 
ai|2,0 >= v^|l,0 > , ai|0,2 >= , ai|l,l >= |0, 1 > 



(35) 



we can determine the realization (on this basis) of the odd generators of osp(l/4; R). 
We obtain 



and 



+ 1 






a = 





A 

dx 





OA 

ax 





(36) 











/ 





\ 




( ° 








1 










X 








Vo 


V2 


j 




[ 


V2x 






(37) 



This imphes (as expected) 

/ \ 

N= 10. (38) 
V 2 / 

The QES BH Hamiltonian is then the sum of the BH Hamiltonian (33) and the QES 
part 

/ 2A((x- 1)^-1) 



Hx^ -2X{x + l) A(^(x-l)i-V^) 

V -V2Xix+l) 

If A = 0, the respective eigenstates and eigenvalues are 

|0,0 >^ E = 0, 
|1,0>±|0,1>^£; = T2, 



(39) 

(40) 
(41) 



|2,0 > -|0,2 



|2, > +|0, 2 > -^{E + 7)|1, l>^E^-^± y^' + 64. (42) 
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If A 7^ 0, the three sectors are mixed and this gives 

XV2{\1, > -|0, 1 >) + (2 - E){\2, > -|0, 2 >) (43) 



for the eigenvalues 



^ = l-|±^\/(7 + 2)2 + 8A2, (44) 



while the states 

ci|0, > +C2(|1, > +|0, 1 >) + C3(|2, > +|0, 2 >) + C4|l, 1 > (45) 

with 

ci = 4:V2{E + 7 - 4)A2 (46) 
C2 = -V2E{E + -f -4:)X (47) 
C3 = -(4^2^ (8-2A2)E- 32) (48) 
C4 = V2{E^ + (7 + '2)E^ + (27 - lOA^)^; - 8A^7 (49) 
are associated with the solutions of 

+ (7 + 2)E^ + (27 - 16 - 12X^)E'^ + (16A^ - 107A^ - 32)E + 128A^ = 0. (50) 

5 Diagonalisation of the Hamiltonian 

In order to achieve this diagonalisation, we will now use a suitable basis of the sub- 
spaces Vn, n = 0, 1, 2. It turns out to be helpfull to write down the matrix elements 
Hij in a simple way for generic values of the number / of oscillators. Along with 
[1, 2], we define the discrete momentum 

k=^ for^=^,^,...,-^ (51) 

The vector containing no quanta is noted 

|0) = [00...0] (52) 

The vectors containing a single quantum are treated by mean of the / vectors of the 
form 

\Mk)) = ^j:{e''Ty-'[lOO...O] (53) 
V / j=i 

and the ^^^^^ vectors containing two quanta are of the form \ijj2,bik)) 

\^2Ak)) =^j:{e^'Ty-'[200...0] 
V/ j=i 

1 ^ 
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(54) 



where the index b takes values b = 1,2,..., and it is understood that there are 
6 — 2 "0" between the two "1" in the different vectors (apart from the case b—1). 

The restriction Hji of H to the invariant vector-space Vq ® Vi ® V2 leads to a 
matrix form 

Hoi 



H 



R 



Hqi Hii H12 







(55) 



where 



iHii)f,^ = -25^^ cos 



5.1 / odd 

Then, if / is odd. 



and 
while 



12^6,1/ = -\/2A5^^5i6 - A (1 + exp — —)5^^5j+ib, 6 = 1,2, 

j=i J 



/ + 1 



22 



-f^22 = -7 - 4 if / = 1 

-7// -%/2g*/^ Oi 
-V2 g// Ow2_3^ -g*// 



O 



-qlf 



-pif J 



\ 

for other values of odd / with 

2i7rp z(/ + l)7a/ ^(/ - l)7r^ /-I /-3 

q — 1+exp — — , p — exp ^ hexp ^ — , i/ — 



(56) 
(57) 

(58) 
/-I 



/ 



/ 



/ 



while // stands for the identity matrix of dimension / and Oip_i^ for the null matrix 
of dimension (|/^ — |/). 

The Hamiltonian further sphts into 1 block of dimension and (/ — 1) blocks of 
dimension 

The different blocks can still be labelled by the momentum k which, along with [1] 
renders the classification of the eigenvalues in function of k possible. Note that this 
contrasts with the pure BH case [1] where there are / blocks of dimension (/ + l)/2. 
The occurence of the supplementary dimensions in the blocks is due to the fact that 
the invariant subspace contains the vectors with zero and one quantum as well. These 
vectors naturally mix with the ones of V2 in the diagonalisation. 
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5.2 / even 

Now, if / is even, 

{Hi2)ijh,u = -\/2X5^Ji^-^/2\5^Ji_^^^-\ (1+exp — ^)5^Jj+ii„ 6 = 1,2, -+1 

(59) 

the second term appearing when // is even only. We also have 

H22^ -7 -4 I if / = 2 (60) 



while 



22 



-7^/ 
-V2g// 



f 

-7 -4 
\ -4 

-V2q*If 
Of 



Olf2_2f 

-q*hf2_2f 



-qlif2_2f 



o 



a/ 



for other values of even / with 



St 



*^(^/2-2/) * (I) 

Of 

2 



(61) 



/ 



2' 2 



+ 1. 



In this last matrix, the notation Om * n stands for the null matrix with m rows and 

n columns while the matrix B is equal to —V2q* Y,j=i ^2j-i,j with em,n standing for 
a matrix where we can find zeroes everywhere except at the intersection of the m*'* 
row and the n*'* column where a 1 is. 

The Hamiltonian further splits into 1 block of dimension | blocks of dimension 
and blocks of dimension 



6 Examples 

We have studied the spectrum of the matrix Hr for a few values of /. As said above 
the different eigenvalues can still be labelled by k and it is possible to plot them on a 
diagram with the momentum k set on the horizontal axis. Drawing such graphics for 
different / and fixed values of A, 7 leads the same pattern as the one of Fig.l of [1]. 
The eigenvalues of are represented on Fig. 1 for A G [0, 0.5] and for 7 = 3, / = 3. 
The label n refers to the number of quanta defined naturally in the A = limit. The 
two lowest lines (labelled n = 2, k = ±1) represent the evolution of the soliton band's 
eigenvalues for the QES-extended BH model. It is clearly seen that these eigenvalues 
stay below the others for a large interval of the new coupling constant A. A similar 
analysis in the cases / = 5, 7 indicates the same phenomenon. 

In order to illustrate the evolution of the algebraic part of the spectrum and the 
splitting of the sohton band in the QES model, we superpose on Fig. 2 the eigenvalues 
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available in the case 7 = 3, / = 7 for two values of A. Here the eigenvalues are 
plotted as functions of i/; the symmetric part (i.e. for v —v) has, of course, to be 
supplemented. The solid-black (resp. dashed-red) hues join eigenvalues corresponding 
to A = (resp. A = 0.5). The various eigenvalues are represented by triangle, square 
and bullet symbols according to the fact that they are related to the n = 0, n = 1 
and n = 2 sectors occuring in the A = limit. The lower curve represents the soliton 
band. The picture clearly suggests the persistence of this band for A > 0. All other 
eigenvalues turn out to be located inside an envelope. To finish, we present some 
detailed calculations of the eigenvalues and of the eigenvectors for / = 1, 2, 3, 4 

6.1 The cases f=l and 3 

If / = 1, the energies satisfy the equation 

+ (7 + 6)^2 + (27 + 8 - QX^)E - IGA^ - 47 A^ = 



which gives for — 3 



A 


El 


E2 


E3 


0.0 


-7.000 


-2.000 


0.000 


0.1 


-7.004 


-2.016 


0.020 


0.2 


-7.016 


-2.061 


0.077 


0.3 


-7.036 


-2.132 


0.168 


0.4 


-7.064 


-2.221 


0.286 


0.5 


-7.101 


-2.323 


0.424 



for the eigenstates 

2\/2A2|0 > -V2XE\1 > +{E^ + 2E - 4X'^)\2 > . 

If / = 3, we have a split : one block of dimension 4 and two others of dimension 3. 
Concerning the first block, the related equation is 

E^ + {^ + A)E^ + (47 - 4 - 18X'^)E^ + (47 - 16 - 12A^ - 167A^)£; + 96A^ - 247 A^ = 
which gives for 7 = 3 



A 


El 


E2 


E3 


E4 


0.0 


-5.372 


-2.000 


0.000 


0.372 


0.1 


-5.389 


-2.043 


0.058 


0.374 


0.2 


-5.439 


-2.159 


0.212 


0.386 


0.3 


-5.524 


-2.314 


0.336 


0.503 


0.4 


-5.645 


-2.484 


0.353 


0.776 


0.5 


-5.801 


-2.649 


0.357 


1.094 



for the eigenstates 

aV6{E + 1)A^|0 > -2V2E(E + l)A|V'i(0) > +(-AE^ + (4A^ - 8)E + 48A^)|V'2,i(0) > 
+V2(E^ + 5E^ + (6 - 14A^)£; - 36A^)|V'2,2(0) > . 
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Figure 1: The energy eigenvalues corresponding to the case / = 3, 7 = 3 are plotted 
as functions of the coupling constant A. 
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Figure 2: The energy eigenvalues corresponding to the case / = 7, 7 = 3 are plotted 
as functions of the momentum v for A = 0.0 and A = 0.5. 
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Concerning the two other blocks, we have the same equation 



+ (7 - 2)E'^ - (3A^ + 27 + + 2 + 7 + 6A^ - 7A^ = 
whose solutions are 

E = 1 , E = -l± y^3(2 + A2) 

or for 7 = 3 



A 


El 


E2 


E3 


0.0 


-3.450 


1.000 


1.450 


0.1 


-3.456 


1.000 


1.456 


0.2 


-3.474 


1.000 


1.474 


0.3 


-3.504 


1.000 


1.504 


0.4 


-3.546 


1.000 


1.546 


0.5 


-3.598 


1.000 


1.598 



The eigenstates are 



> -^|^2,i(±f ) > +^(l±^V^)|V'2,2(±f ) > 

for E — 1 and 

27r 271 271 

-3(lTi\/3)A|V'i(±y) > +\/2(1tzv^)(E-2)|V'2,i(±y) > +2(^ + 1)|^2,2(± y) > 

for the other energies. As is clear from above, we observe a twofold degeneracy : the 
degenerated states are complex conjugated. 

6.2 The cases f=2 and 4 

Some examples : 

If / = 2, we have one block of dimension 4 and another one of dimension 2. The 
energies of the first block satisfy the equation 

+ (7 + 2)E^ + (27 - 16 - 12X'^)E^ + (-32 + IGX"^ - 10^X'^)E + 128A=^ = 



which gives for 7 = 3 



A 


El 


E2 


E3 


E4 


0.0 


-5.772 


-2.000 


0.000 


2.772 


0.1 


-5.782 


-2.029 


0.039 


2.772 


0.2 


-5.813 


-2.110 


0.150 


2.773 


0.3 


-5.865 


-2.227 


0.318 


2.775 


0.4 


-5.939 


-2.364 


0.525 


2.777 


0.5 


-6.034 


-2.507 


0.761 


2.780 
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for the eigenstates 

A{E - 1)A2|0 > -V2E{E - l)A|V'i(0) > -{4E^ + (8 - 2X'^)E - 32A2)|V'2,i(0) > 
+{E^ + 5E^ + (6 - lOA^)^; - 24A2)|V'2,2(0) > 

The energies of the second block satisfy the equation 

E^+{-f- 2)E - 27 - 2A2 = 



which gives for 7 = 3 



A 


El 


E2 


0.0 


-3.000 


2.000 


0.1 


-3.004 


2.004 


0.2 


-3.016 


2.016 


0.3 


-3.036 


2.036 


0.4 


-3.063 


2.063 


0.5 


-3.098 


2.098 



for the eigenstates 

V2A|V^i(7r) > +(2 - ^)|V2,i(7r) > . 

If / = 4, we still have a split : one block of dimension 5, another one of dimension 
4 and two of dimension 3. Concerning the first block, the related equation is 

E^ + (-f + 2)E'^ + (27 - 16 - 24A^)£;^ + (-87 - 32 + ^2\^ - 22-fX'^)E^ 
+(-167 + 256A^ + 167A^)£; + 1287A^ = 



which gives for 7 = 3 



A 


El 


E2 


E3 


E, 


E, 


0.0 


-5.191 


-2.000 


-1.317 


0.000 


3.509 


0.1 


-5.214 


-2.066 


-1.307 


0.078 


3.509 


0.2 


-5.282 


-2.228 


-1.288 


0.289 


3.509 


0.3 


-5.398 


-2.429 


-1.273 


0.590 


3.510 


0.4 


-5.563 


-2.631 


-1.263 


0.947 


3.510 


0.5 


-5.778 


-2.814 


-1.257 


1.338 


3.511 



for the eigenstates 

4y/2{E -4:){E + 3)X'^\0 > -V2E{E - 4){E + 3)A|^i(0) > 

+ {-8E{E + 2) - 2(-64 + E{E - 8))A^)|V'2,i(0) > +2\/2(E + 3){-E{E + 2) 

+(£; + i6)a2) 1^2,2(0) > 

+(128A' + E{E + 2)(-8 + E{E + 3) - 22A2)) 1^2,3(0) > . 
Concerning the block of dimension 4, we have 

E(E^ + (7 - 2)E^ - (4A^ + 27)^; - 27A^) = 
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which gives for 7 — 3 



A 


El 


E2 


E3 


^4 


0.0 


-3.000 


0.000 


0.000 


2.000 


0.1 


-3.004 


-0.010 


0.000 


2.014 


0.2 


-3.016 


-0.039 


0.000 


2.055 


0.3 


-3.036 


-0.084 


0.000 


2.120 


0.4 


-3.065 


-0.142 


0.000 


2.206 


0.5 


-3.101 


-0.209 


0.000 


2.311 



The eigenstate is 

IV'2,2(7r) > 

for the null energy and the other states are 

V2{E + 3)X\x/ji{7r) > -2A2|^2,i(7r) > -{E^ + E - 2X'' - 6)\ij2M > 

in what concerns the three other energies. 

Concerning the two final blocks, we have the same equation 



and for 7 = 3 



A 


El 


E2 


E3 


0.0 


-4.000 


0.000 


1.000 


0.1 


-4.009 


0.005 


1.004 


0.2 


-4.036 


0.020 


1.016 


0.3 


-4.080 


0.043 


1.037 


0.4 


-4.140 


0.072 


1.067 


0.5 


-4.215 


0.107 


1.107 



The eigenstates are 



(l + E)A|^i(±-) > -x/2(A2-i5)|^2,i(±-) > --{l±z){E' + 3E-2X')\i^2,2{±^) > ■ 

As is once again clear from above, we observe a twofold degeneracy : the degenerated 
states are complex conjugated. 

7 Concluding remarks 

The exactly solvable models of the BH type can naturally be generalized to a family 
of /-body quasi exactly solvable, translation invariant Hamiltonians. These general- 
isations depend on one (or more) new coupling constant (s). Here we have put the 
emphasis on a peculiar property of the spectrum in the generalized models : when 
plotted with respect to the discrete momentum k, a line of eigenvalues (one for each 
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values of k) appears splitted from the rest of the spectrum, forming a sohton band. 
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